VIBRATIONS OF A BEAM BETWEEN STOPS: CONVERGENCE 
OF A FULLY DISCRETIZED APPROXIMATION 
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Abstract. We consider the dynamics of an elastic beam which is clamped 
at its left end to a vibrating support and which can move freely at its right 
end between two rigid obstacles (the stops). We model the contact with Sig- 
norini's complementary conditions between the displacement and the shear 
stress. For this infinite dimensional contact problem, we propose a family of 
fully discretized approximations and their convergence is proved. Moreover 
some examples of implementation are presented. 
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1. Description of the problem 

We consider a beam which is clamped at its left end to a vibrating support and 
which can move freely between two rigid obstacles at its right end (see figure 1). 

The longitudinal axis of the beam coincide with the interval [0, L] and we denote 
by u{x,t), {x,t) G (0,L) X (0,T) the vertical displacement of a point x belonging 
to this axis. We assume that the material is elastic and the motion is planar. We 
denote by a the shear stress given by 

d'{x,t) = —k Uxxxj k = — ■ 

pS 

where p and E are the density and the Young's modulus of the material and S 
and / are respectively the surface and the inertial momentum of the section of the 
beam. Then, under the assumption of small displacements, the motion is described 
by the following partial differential equation 

Utt -^x = f 

where / is the density of external forces. 




Figure 1. The physical setting 
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The beam is clamped at its left end so 

M(O,i)=0(t), i2,(0,t) = 

where (p describes the motion of the vibrating support. At its right end the beam 
can move freely between two obstacles, called "the stops", so we have 

9i < u{L, t) < g2, Uxx{L, i) = 

and we assume that 51 < < 32- When the beam hits one of the two stops, the 
stress is in the opposite direction of the displacement and we obtain the following 
Signorini's conditions 

a{L,t)>0 \iu{L,t)=gi, 
alL,t)<0 iiu{L,t)=g2, 
a{L,t) = if gi < M(L,t) < 52- 

These relations can be rewritten as follows 

~d-{L,t) e ^^^:[g^^g^]{u{L,t)) 

where '>Plgi,g2] the indicator function of the interval [(71,32] and dtp^g-^ g^j is its 
subdifferential ([101). 

In order to deal with homogeneous boundary conditions at a; = 0, we consider a 
new unknown function u defined by 

u{x,t) = u{x,t) — h{x)(j){t), 

with 

^(^)-i-Kz) +3(z) -3(1) • 

The mechanical problem is now described by the system 

utt + k^Uxxxx = f in (0, L) x (0, T) 
uiO,-)^UxiO,-)^Uxx{L,-) = in (0,T) 
u{L,-) e [91,92], Uxxx{L,-) e d^[g^^g^]{u{L,-)) in(0,r) 

with f{x,t) = f{x,t) - h{x)4>" (t) - Ph(-^^x)(j)it) for aU {x,t) e (0, L) x (0,T). We 
complete the model with the initial conditions 

u(-,0) = uo, ut{-,0)=vo in(0,i). 

As usual in mechanical problems with unilateral constraints we cannot expect 
classical solutions since the velocities may be discontinuous. So we look for weak 
solutions. For this purpose we introduce the following functional spaces 

H = L2(0,L), V ^ {we H^{0,L);w{0) = Wx{0) ^ O}, 

H = {w e L^iO,T-V);wt e L^{0,T;H)], 
and the convex set 

K ^ {w£V;9i<w{L)<g2]. 

We denote by (., .) and |.| the canonical scalar product and norm of H . Let a be 
the following bilinear form 

a{u,v) = I k'^UxxVxxdx V(?i,w)Gy^. 
-'0 
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We may observe that a defines a scalar product on V and the associated norm, 
denoted is equivalent to the canonical norm of H^{0,L) on V. The weak 

formulation of the problem is then given by the following variational inequality 



Vw e nnL^{0,T;K) such that w{-,T) = m(-,T). 

For this problem an existence result has been obtained by K.Kuttler and M.Shillor 
by using a penalty method. 

Theorem 1.1. (15] j Assume that f e L^{0,T;H), uq e K , vq e H. Then there 
exists u eTIC] L^{0, T; K) such that problem (P) is satisfied and u(-, 0) = uq. 

It should be noted that, as far as we know, uniqueness remains an open question. 

For the computation of approximate solutions, the penalty method which is 
introduced as a theoretical tool to obtain existence in [S] could appear as an inter- 
esting technique: the Signorini's conditions are replaced by a normal compliance 
law 

a{L, t) = — i [max(u(L, t) — 52, O) — max((7i — i), 0)1 , - >> 1 

£ £ 

which leads to a system of partial differential equations depending on the penalty 
parameter e. From the mechanical point of view 1/e. can be interpreted as the 
stiffness of the stops which are not assumed to be perfectly rigid anymore. From 
a numerical point of view, for large values of 1/e, we have to solve a stiff problem 
which is expensive (|2I)- Moreover the dynamics of the system may be complex 
(see for a periodic forcing) and the approximate motion could be quite sensitive 
with respect to the value of 1/e. (see |S| for an example in the case of a simplified 
model of vibrations, see also ^). 

In order to avoid these difficulties, we propose to deal directly with the unilateral 
boundary condition by solving a complete discretization, in both time and space, 
of the variational inequality (P). 

From now on we will consider the more general case of a convex set K given by 



where 51, gi are two mappings from [0, L] to R such that there exists <? > such 



Problem (P) then describes the vibrations of an elastic beam between two longitu- 
dinal rigid obstacles. 

The paper is organized as follows: in the next section we introduce the fully dis- 
cretized approximation of the problem, then in section|21 we prove its stability and 
convergence and finally, in section 01 we present some examples of implementation. 

Let us observe that the convergence result yields also an existence result for the 
more general case that we consider here. 




K={w(^ V;gx{x) < w{x) < g2ix) Vx e [0,L]} 



that 



gi{x) <-g<0<g< g2{x) Vx- G [0, L]. 
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2. Discretization 

Let us assume now that the assumptions of theorem ll.il hold i.e / G i^(0, T; H), 
uq £ K and wq S H. 

For all h G M*^ we consider a finite dimensional subspace of V such that, for 
all V Cz there exists a sequence {vh)h>Q such that 

- v\\v ->/i-.o 0, w/i eVh yh> 0, 

and we denote by Qh the projection onto Vh respectively to the scalar product 
defined by a on V. The compact embedding of V into H^{0, L) implies that there 
exists a sequence {•^h)h>Q such that 

For all /i > we define Kh^ Kf] Vh- 

Let A'' e N* and Ai = T/A^. We propose the following family of discretizations 
of problem (P): For all n e {1, . . . , - 1}, find ul+^ G Kh such that 



^ >\ +a{pul+^ + {l-2p)ul+pul-\ii,h-ul+^) 

> {f3r+^ + (1 - 2/3)/" + pr-\wh - ul+^) Vwh e if/. 



with 



1 



At 



(n+l)At 



where /3 is a parameter belonging to [0, 1/2]. We choose w° and uj^ in such that 



|u^||y)^^Q remains bounded and 



(2.1) lim ||<-uo|k + 



Ut -u 



h 



At 



0. 



We may observe that problem (P^/^) (/? G [0, 1/2]) can be rewritten as 
f Find G Kh such that 

\ anp{ul+\wh - ul+') > Lnp{wh - ul+') ywh G Kh 

with 

anpiuh.Vh) = {uh,Vh) + At^f3a{uh,Vh), 
L^0{vh) = Ai2(/3/«+i + (1 _ 2/3)/" + /3r-i,^;„) + (2< - 
-Ai2a((l-2/3X+/3<-\f„) 

for aU {uh,Vh) G 

By an immediate induction on n, we obtain that L„/3 is linear and continuous 
on Vh, and it is obvious that a„^ is bilinear, symmetric, continuous and coercive 
on Vh- Thus the existence and uniqueness of u^^^ follows. 
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Remark 2.1. This family of discretizations is inspired by Newmark's algorithms 
of parameters 7 = 1/2, /? e [0,1/2] ( 4 ). Indeed, li K = V (i.e gi{x) = -00, 
52(2;) = +00 for all X G [0, L]), the mechanical problem is described by the following 
system of partial differential equations 

utt + k'^Uxxxx = f in (0, L) X (0,r) 
u(0, •) = Ux{0, ■) = Uxx{L, •) = Uxxx{L, •) = in (0, T) 

with the initial data {uo,vo), and {Pl^g^) reduces to 



find e Kh such that 



<+i-2< + <-i 



u.;,^ + a {Pul+^ + (1 ~ 2/3X + Pul-\wh) 



A<2 

(/3/"+i + (1 - 2/3)/" + Pr-\wh) Vwh e Vh 

which is simply a Newmark's scheme of parameters 7 = 1/2, /3 for the previous 
system. 

3. Convergence 

Since the proposed discretizations are inspired by Newmark's methods which 
stability depends on the value of /3, we may expect the same kind of result for 
i-^hp^)- More precisely, for P g [0, 1/2) we obtain the following conditional stability 
property: 



Proposition 3.1. Let f3 E [0, 1/2), h > and Kh be defined by 

a{uh,Uh) 



sup 



Let a e (0,1) and Nh G N* be such that 

T ' ' 



(3.1) 



< min 2. 



1 - a 



Kh{\-2f3) 



, a 



Then there exists a constant depending only on the data, C{f,Uo,Vo), such that for 

T 

aUh>0 and for all N > Nh (i.e At < Ath = — / 

Nh 



a 



,,"+1 _ ,,n 



At 



+ /3a«, + Pa{ul+\ul+') < C{f, uo,vo) 



for all n € {1, . . . , TV - 1}, where ) i<n<w-i are the solutions of problems 

{Php^)l<n<N-l- 

Remark 3.2. An estimate of Kh in the case of a P3 finite element space discretization 
is given in the Appendix. 

Proof. Let rt e {1, . . . , TV — 1} and choose Wh — u^^^ as a test-function in [Ph^^)- 
We get 



n+i _ r, n , n-i 

a^ u 



+a(/3<+i + (l-2/3)< <-!-<+!) > 
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where 

= + 2/3)/" + /3r-i. 

The two first terms can be rewritten as follows: 



Ai2 



)- 




2 






At 




At 



and 

a (/3<+i + (1 - 2/3)< + f3ul-\ul-' - ul+') = (1 - 2(3)a{ul~\ul) 
+(3a{ul-\ul-') (1 -2/3)a«,<+i) -/3aK+\ 

Hence, for all n e {1, . . . , iV — 1}, we have 



,n+l „,n 



At 



+ (1 - 2/3)a«,<+i) +/3aK+S<'+i) 



< 



At 



+ (1 - 2/3)a«-\ <) + Pa{ul~\ul-') + {g 



and with a discrete integration 



."+1 o," 



< 



At 



At 



+ (1 - 2f3)a{ul, + /3a«+\ + /3a«, <) 

+ (1 - 2/?)a«, <) + /3a«, ul) + /Ja^, 7.°) 



0=1 



Using the same techniques as in [TJ, we define 



R{uh,Vh) = (1 - 2(3)a{uh,Vh) + 



Uh - Vh 



At 



We observe that 

1 - 2/3 



R{uh,Vh) 
1 - 2/5 



^ 1-2/3 , 



> 



a{uh + Vh,Uh + Vh)+ [l- KhAt 



4 

,1-2(3 



Uh - Vh 



At 



Uh - Vh 



At 



and with assumption H3.1() . we infer that 

N 1-2/3 , 
^(w/i, Vh) > — - — a{uh + Vh, Uh + Vh) + a 

It follows that 
(3.2) 



Uh - Vh 



At 



y{uh,vh)€V^. 



n+l 



At 



< R{ul ui) + Pa{ului) + /3a«, <) + J^ls^^At + ^ 

p— 1 p—0 



p+1 p 



At 



At. 
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Since a — At > a — Ath > 0, Gronwall's lemma implies that 



E 

with 



At 



< 



At 



exp 



nAt 
t - At 



E 



kp exp 



{n-p)At\ 
a -At J 



p 

R{ul ul) + Pa{ul ui) + Pa{ul ul) + V j^T Ai 



-/c|2, 



At 



Since / € L^{0, T; H), we infer that the right hand side of (|3.2|l remains bounded 
by a constant which depends only on the data (/, uq, vq). □ 

We may observe that the lack of stability is due to the terms (1 — 2f3)a{u]l, ""/J^^) 
and (1 — 2/3)a(u5J~^, wJJ). For the case /3 = 1/2, this difficulty does not occur and 
we obtain an unconditional stability result: 



Proposition 3.3. Let (3 — 1/2. Then there exists a constant depending only on 
the data, C{f, uq, vq), such that for all h > and for all N > 1 
, 2 



n+l 



At 



+ ^a{ut ul) + iaK+\ < Cif, no, vo) 

,n+l 



for all n £ {1,...,A^— 1}, where (wJJ )i<n<Af-i are the solutions of problems 

iPhf3^)T-<n<N-l- 



We define now an approximate solution wj^pj (/9 G [0, 1/2]) of problem (P) by a 

"+1 of fP"+ 



linear interpolation of the solutions u^'^^ of (Phg^)- More precisely, for all h > 



and iV > 1 



t)^^n (.n+l)At~t +, t-nAt 



for aU t e [nAt, [n + l)At\ , < n < iV - 1. 

Let a e (0,1) and Nh be defined by condition l|3.1|l if /3 e [0,1/2), otherwise 
let Nh — 1 for all h > Q. The previous stability results imply that there exists a 
subsequence, still denoted {uf^ j^)h>o,N>NhJ and u £ W — {w e L°"{0,T;V),Wt G 
X°°(0,r;iJ)} such that 



u 





h,N 



weakly* in L°°{Q,T;V), 



^ ^ weakly* in L°°(0, T; H). 
dt dt V 1 . ; 

With Simon's lemma ('|12j) we infer that W is compactly embedded in C°([0, T];H'^{0, L)) 
and we know also that W C C°'i/2([o,L] x [0,r]) (ill|). It follows that, possibly 
extracting another subsequence, we have 

uI j, ^ u strongly in ^^([O, T]-H\Q, L)), 
and thus u belongs to i^(0, T; K) and 0) = uq. 



Let us prove now that u is a solution of problem (P). 
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Theorem 3.4. Let (3 E [0,1/2] and let Nh be defined by condition iS. 1]) if f3 ^ 
1/2, otherwise Nh — 1 for all h > 0. The sequence of approximate solutions 
(''^h N)h>o,N>Nh admits a subsequence which converges weakly* in W to a solution 
of problem (P). 



Proof. We consider now the converging subsequence of {uf^ j^)h>Q,N>Nh' still de- 
noted iu1^j^)h>o,N>Nn - Let w e nn L'^{0,T;K) such that ?D(-,T) ^ u{-,T). We 
will prove that 

{ut{-,t),wt{-,t) - ut{-,t)) dt + / a{u{-,t),w{-,t) - u{-,t)) dt 

T " 
>{vo,w{-,Q)-uo)+ I {f{-,t),w{-,t)-u{-,t)). 
Jo 

As a first step we have to construct a well-suited test-function . 
Let e e (0, T/2) and (j) he a. C°°-function such that 

" < (f>it) < 1 yte [o,r], 

0(i)=O yt e[T -3e/2,T], = 1 Wt £ [0,T - 2e]. 

We denote w — {1 — <j))u + (pw. Since K is convex, we have immediately w S 
nnL^{0,T;K) and w{-,t) = u{-,t) for all t G [T-3£/2,T]. 

Let 77 G (0,e/2) and ^ G (0, 1). Following the same ideas as in jT^ we define 

w'i't' by 

(3.3) w''-''{-,t) ^u{-,t) + - J {{1- fi)w{-,s)-u{-,s))ds Vt G [0,T-£/2]. 

Since u d W and w G 7Y, we have immediately w''''^ — u G C"([0, r];F), wj*'^ G 
L^{0,T;H) andu;''^^ G L°° (O, T; nC° ([0, T]; iJi(0, L)) . Moreover we can choose 
77 such that w^'^ satisfies strictly the constraint. More precisely, for all t £ [0, T — 
e/2] and for all x E [0,L] we have 

I rt+n Y rt+ri 

w^''^(x,t) = — / (1 — fi)w(x, s) ds + u(x,t) / u(x,s)ds. 

V Jt V Jt 

The first term of the right hand side belongs to [(1 — //)gi(a;), {1 — ii)g2{x)\ with the 
convention that (1 — ii)gi(x) = gi{x) [i — 1, 2) if gi{x) G {+00, —00}, and recalling 
that u G C"^i/2([0,L] X [0,T]) we have 



1 



— ~ [ Wix, t) — u{x, s)\ds< — 



nix, t) / u{x, s) ds 

V Jt V Jt 

where Cq is the Holder continuity coefficient of u. 
Thus, choosing 77 such that 

ensures that 

(3.5) .91(2;) + |.g < w^^^ix,t) < 32 (x) - ^g 

for all i G [0, T — e/2] and for all x G [0, L], with the convention that gi{x) ± ^g — 
gi{x) {i = 1, 2) if gi{x) G {-f^oo, -00}. 
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Now, we assume that At < ^ and, for all n £ {1, . . . , — 1}, we define the 
test-function by 

^■^■^> - \ if nAt>T- e. 

We have to check that belongs to Kh- 

Lemma 3.5. There exist hi > and N'f^ > Nh such that, for all h G (0,/ii) and 
for all N > N'l^, we have 

wl&Kh Vjie {!,..., TV- 1}. 

Proof. First of all it is clear that ui^' € Vh and G Kh if nAt > T — e. Otherwise, 
when nAt <T — e, we rewrite as follows: 

, -^ < = <jv(-' (" + 1)^*) - + + + 1)^0 - <-,nAt) 

^ ■ ' +w■^■^'{■,nAt) + {Qh - I){wi'^'{-,nAt)-u{-,nAt)). 

We already know that j^)h>Q.N>Nh converges to u strongly in C°([0, T];H^{0, L)) 
and u e C0'i/2([o,i] x [o',r]), thus 



(3.8) sup \u{x,{n + l)At) -u{x,nAt)\ <CoVAi 

xe[o,L] 

and 

(3-9) ^snpJul^{x,{n+l)At) ^u{x,{n + l)At)\ < Ci\\ul^-u\^^ 

where Ci is the norm of the canonical injection of H^{0, L) into C''([0, L]). More- 
over 
(3.10) 
sup. 



■0 L] I (Q/. - I) (w"''' (a^, «Ai) - u{x, nAt)) \ 
{Qh - I){w'>'^{-,nAt) - u{-,nAt)) \\^,^^^^^ < Cij^w^''^ - u\\^^^^^^,^y 



By choosing hi and N'f^ > Nh such that 

|g > CoVAi + Ci\\ul^ - "llco([o,T];ffi(o.L)) + C'n/.llu;"^^ - u||Lo=(o,T;y) 

for aU h e (0, /ii) and At = ^ with iV > iV,;, relations and (|T7|l - l|TTn|) imply 
that Wh{x) belongs to [gi{x), 52(2;)] for all x e [0, L] if nAt < T—e which concludes 
the proof. □ 



Let us choose now Wh = w"^ in (Ph/^^), 1 < rt < iV — 1 and define N' = 
With a discrete integration we obtain 



T-e 
At 



At 



N' 



N' 

(3.11) < J2 <(^<^' + (1 - 2/5)< + < - At 

ri=l 



n=l 
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and we have to pass to the hmit in each term as h and At tend to zero. RecaUing 
()2.1(l we immediatly infer that 

(^^^' - <) - (-0, W^-^^{;0) - U(, 0)) . 

Then, we rewrite the second term as fohows 

N' N' 



n— 1 n— 1 

But 



Af' N' »(„+i)At 
„_i 1 JnAi 



/.(n+l)At 

E 

n=l 



nAt 



(/(•, s), K''' - u){;nAt) - - s)) ds 



n=l 

Observing that 

\{Qh - - u){;nAt)\ < - u|U-(o,T;y) 

for aU n £ {1, . . . , N'}, we obtain 



AT' 



(^^'^ - ^)((^''^'' - u)i;nAt)))At 



N' 



<J2lk\fV'' -u\\L^iO,T;V)At 



< VT'jh\\f\\mo,T:H)\\w'^''' - u\\l'=°{o,t-v) 0. 

Moreover, with the definition of w^'^, we have 
(3.12) 

_ u){-,nAt) - (w''''^ - u)(-,s))||^ 



1 

< - 



V J nAt 



{{I- fi)w- u) {-,(7)11^(1(7+ / ~ fi)w - u){-,a)\\yda 



2y/ s — nAt 
< 11(1 - WW - u||L2(o,T;y) 



riAt+r; 



for aU n e {0, • • • , N'} and s e [nAt, (n + l)At]. 

If we denote by C the norm of the canonical injection of (V, \\ ■ \\v) into | • |), 



we get 



N' ^(n+l)At 



E 



(/(•, s), {w^^'' - u){; nAt) - [w^^^ -u)i-, s)) ds 



< VAt 11(1 - - 'u|lL2(o,T;y)||/||L2(0,T;H) ^ 0. 
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Finally, since / E L^{0,T;H), w''■^'-ue C°{[Q,T];V) C L^{0,L;H) and |T-e- 
{N' + l)At\ < At, we may conclude that 



N' ^(n+l)At 
lim V / {f{;s),{w^''^~u){;s)) 



ds 



n=l 
T-e 



nAt 



and 



N' 



T-e 



Moreover, relation (|3.12() implies that 



2VAt| 
77 



1(1 



(3.13) ||(w"-^-u")-(w"-u"+^)||^ < 
for aU 71 e {1, . . . , iV' + 1}. It follows that 

AT' + l 

^ (/" - /"-I, (w"-i - u") - (lo" - At 

n=l 

Finally we observe that 

\Wh - <+'lk < ll^""^^ - "IU-(o,T;y) Vn e {0, . . . , N'}. 

Hence 

|(/° - f\w' ~ u^)At\ < 2CVAt||7i;"^^ - "IU==(o,T;y) II/IIl=(o,T;//) 
and we may conclude that 

N' T-e 

^(5",<'-<+i)At^ / (/(,5),(7«"-'^-K)(.,5))ds. 

n=l 



V) 



Let us study now the convergence of the first term of the right hand side of 
l)3.11|l . First we rewrite it as follows: 



N' 



a\^-^ul + Pulw1~ul]At 



(3.14) 



- E « u^^r' + -^<' w - "^') - K"' - ) At 

n=l ^ 
JV' 



E 

n=l 
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With the propositions 13.11 and VA.'M wc know that (ll^f jvIU°°(o,T;\/)),j>q ^>jv^ 
bounded independently of h and At, thus 

1 



IS 



< §ll"h,Ar||L-(0,T;y)||w''''" - u||l~(0,T;V') 



and, with l|TT^ 



< ll<,Jvlli°°(O.T;V)||(l - - w| 



77 



L2(0,T;y)- 



FinaUy, we rewrite the last term of H3.14|l as follows 

N' / „ + i , „ N' ^(„+i)At 



E 



n+l I n 



WAt 

a{ulj^{-, s), Qh{{w'''^ - u){;nAt)) - Qh{{w'^'^ - 5)) ds 

nAt 







E 



n=l 
rT-e 



Jo 

The two first terms can be estimated by 
and, with estimate (|3.12|) we have 

^(ri+l)At 



2TVKt. 



< 



V 



l!"/l,A'IU~(0,T;y)||(l - ^^)W - M|L2(o,T;y)- 



Finally, recalling that Qh{w'^''^ — u) is the orthogonal projection of w''"'^ — u on 
Vh respectively to the scalar product defined by a on we obtain that 



T-E 



Jo 

and the weak convergence of Uh to u in L^{Q, T; V) allows us to conclude that 

N' .T-e 

^ a(/3<+i + (1 -2/3K+/3<, At ^ / a(u(., s), (u;"'^ - s)) 

„=i -^0 



ds. 



There remains now to study the convergence of the last term of (13.11(1 i.e 



(3.15) 



E 



At 



At 



^ At. 



In order to simplify the notations, let us define 

(w'''^ - t + At) - (w''^^ - u)(-, t) 



'4'Ati-,t) = 



At 



, vte[o,r-e]. 
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We rewrite (|l-i.l5|l as follows: 



E 

n=l 



At 



At 



+ E ^ A, ^ (Q/^ - ^)(V^A*(-, - l)At)) At 

n=l ^ ^ 
„=l"'(«-l)At 



At 



,nAt ^^n n-1 



E, 



V'At(-,(n-l)Ai)-i^At(-,i)) rf< 

"""" ' ^At dt. 



The first term, which can be interpreted as a boundary term aX t — T for the 
discrete time integration, can be estimated by 



" " " ,,,JV' _ ,,^' + 1 



Ai 

C 

< — max 

77 l<n<Af 



At 



N'At+ri 



N'At 



^J.)w{■,s) -u{-,s)\\ ds. 



But iV'At > T - 3e/2, thus «;(•, s) = s) for aU s e [iV'At, iV'At + r]] and 



N'At+T] 



N'At 



N'At+r] 



fi)w -u){-,s)\\yds = / //||u(-,s)||^ds < /i77||u||i^(o,T;y)- 



TV' At 



Since proDOsitions 13 . ll and 13 . 3l imolv that max 



y 

n-l 



At 



is bounded indepen- 



dently of h and At, we infer that there exists a constant C such that 



,N' + 1 „,N' 



At 



N' N' + l 



< CV||w||l~(o,T;V)- 



For the second term we perform the same kind of computation: 



N' / „ n-l 



E 



n=l 
N' 

n=l 



At 



^,iQh-I){^At{;in^l)At))]At 



At 



liw"''' - u){-,nAt) - (w"'^ - u){-, (n - l)At)| 



< — max 

77 l<n<Ar 
nAt+rj 



At 



N' 

E 

n=l 



lAt 



(n-l)At 



|((1 - fi)w ~ u){-,s)\\^ds 



+ 



{n-l)At+ri 



|((1 - fi)w -u){-,s)\\yds 



< 2 max 

77 l<n<JV 



At 



|(l-A*)w^-"|L2(o,T;y)- 
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Recalling that {"fh)h>o converges to zero, we obtain 

E ( ' (^'^ - ^)(V'A*(-, (r^ - l)At)) ) At - 0. 



In order to estimate the third term, we transform ipAt{-, {n — l)Af) — tpAt{-,t) as 
follows: 



^Af(-,(n-l)At) -7/.Ai(-,t) 



^ [C-DaM^^ - - t^) (•, 5 + At) - ((1 - - (•, s)) 



ds 



(n-l)At+?; 



((1 - - m) (•, s + At) - ((1 - fj.)w -u){; s)) 



p(n- 

+ 

Jt+ri 

1 rs+l^t 

= — / / ((l-/i)wt -Ut)(-,cT)do-rfs 

riL^t 7(«-i)At 



ds 



(n-l)At 

(n-l)At+?) ^s+At 



+ 



Hence 



((1 - fx)wt - ut){-,a)dads. 



|V'At(-,(n- l)At) -VAt(-,i)| < ^^^-^^^7=J^^||(l-/u)«;t-Ut||^2(Q_T;_H') 



?7VAt 



and 



N' ,nAt /yn_yn-l 



E 



n=l 
N' 

<E 



(n-l)At 
At2 



^ ,^At(-,(n-l)At) -^At(-,^)) rfi 



At 



1(1 - IJ,)Wt - Ut\ 



L^{0,T-H) 



^ TV At 

< max 



r] l<n<N 

Finally, we observe that 



At 



|(l-/x)«;t-Wt||i2(o,T;H) 



N' ^„At 
n=i-'(n-l)At 



At 



T-s ( „N' 



-f 

Jo 



dt 



i;t),i;At{;t) dt 



JN'At \ 



At 



-,V'At(-,i) dt. 



Since {1~ n)w-ue L^{0, T; V), and 



I / H+ri+M pt+At \ 

V'At(-,t) = ^ ( y {{1- fj,)w-u){-,s)ds- J {{1- tJ,)w-u){-,s)dsj 

((1 - ll)w - u){; t + 77) - ((1 - ll)w - u) {; t) 



we obtain that 

V'At(-,i) -^At^o 
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Strongly in (0, T - e;V). Since - — — ^ — converges weakly to — in i^(0, T; V), it 

at at 

follows that 



{■■.t),iJAti-,t) dt 



\ dt J 

'du ((1 - ^l)w - U) {; ^ + 7?) - ((1 - fi)w - U) {; t) \ 

d^'' 77 j 

Moreover, for alH e [0, T - e] 

\i^At{-,t)\<^\\{w^-^^^u){-,t + M)-{w'^-^^-u){-,t)\\^ 
< — ^=1 (I ~ ii)w — u\\ 



?]VAt 



\L^{0,T;V)- 



Thus 



< 



[ ' .M;t)) dt 

N'At y At J 

T] ^' ^^L^iO'^y) l<n<N 



At 



Taking into account the previous convergence results, we obtain 

\f{-,t), K-'^ - U){; t)) dt + {Vo, K-'^ - U){; 0)) 



< 



a{u{-,t),{w'^''' -u){-,t))dt 
du^^ ((1 - ^)n^-u)i-,t + rj)- ((1 ~ ii)w - u) (■, t) ^ 



dt - 77 

+ CVII'«||L-(0,T;y) 

for aU e e (0, T/2), ^ e (0, 1) and 77 e (0, e/2) satisfying 

Rewriting ((1 — /z)7i; — w) (•, < + 7/) — ((1 — /i)?xi — ?/) (•, t) as 

rt+n 

((1 - /x)7«f - ut){-,a) da 



and recalling that (1 — fi)wt ~ ut E L^{0, T; H), we obtain that 

- / ((1 - fi)wt - Ut) {a, •) da ^,,^0 (1 /^)wt - Ut 

strongly in (0, T^e;H). Similarly, since {l~^j.)w~u e 1^(0, T; V)nC"{[0, T] ; H) 
we have 

1 /■*+'' 

strongly in (0, T — e;V), and 

(w''^^ -u)(-,0) ^,,^0 (1 -/i)w(-,0) - u(-,0) strongly in i7. 
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Thus, when rj tends to zero, we get 

rT-e 



[ {f{;t), ((1 - ^i)w - U){;t)) dt + [Vo, ((1 - fi)w -u){; 0)) 
j-T-e 

w — u){-, t)) dt 

r^-^ fdu \ 
- J { ~ ~ ) '^^ + C''/^ll"IU~(o,T;y)- 

Finally, we can pass to the limit when and e tend to zero and, observing that 
w — u = <f){w — u), we may conclude the proof. □ 

4. Finite element implementation in (P^^"^) 

We present now some simulations when the contact with the stops takes place 
only at the right end of the beam, i.e 

gi (x) = —00, g2 {x) = +00 Vx G [0,L), 
91 {L) = -g, 92 {L) = g, 

with g a positive real number. We use the well-known Hermite piecewise cubics 
as basis functions for the space discretization. So we consider a partition of the 
interval [0,L] into J subintervals of length h, i.e. xq = 0, Xi = ih, xj = L. At 
each node Xi, we associate two Hermite piecewise cubics (p2i-i and (fi2i {i = I, J) 
defined by 

V2i-i G P3, V2i-i (xj) ^ Sij and (/?2i-i (xj) = for 1 < j < J. 
V2i e P3, V2i (xj) = Sij and ifi2i (xj) = for 1 < j < J. 

Thus, if 

J ,/ 

W = ^ U2i-liP2i-l + ^ U2ilf2i, 
1=1 i=l 

the coefficient W2i-i gives the value of u at node Xi and the coefficient U2i gives 
the value of the derivative of u at node Xi , 1 < i < J . Hence, for any u & V ^ the 
interpolate Uh is given by 

J J 

i=l i=l 

We consider the following finite dimensional subspace 

Vh = span{v?i, </92, • • • , ^2j-i,^2j} C V. 
Let AT > 1, /? e [0,1/2] and «+^)o< „<jv-i be the solutions of problems 

hf3 

2.7 



(P^^^)o<n<JV-i- Since belongs to Kh = Vh D K we have 



= E^rVi with ul+\L)=u-+\ e [-5,5] 
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can be rewritten as follows: 



. — Tn>2J-2 



■n+l 



[^5> 5] X K. such that, for all w E Kh 



Ai2 



~i7l — —n-\-l\ 



where 



+ (1- 2/?) /"+/3/"-\(p,) I 



1, 



.,2J 



and and S are respectively the global mass and stiffness matrices. The previous 
inequality is also equivalent to the differential inclusion: 



M 



TJl + l 



2u" 



Ai2 

3 G" 



+ 5 + (1 - 2/3) u" + 



which can be rewritten as 

(4.1) {M + AfpS) + At^dt^K^ (ll"+i) 

with i^" = (2A^ - At^{l - 2/3)5) u" - {M + At^pS) u" 



3 

More precisely, if we consider a single beam element [xi,a;i+i] of length h (i 
0, . . . , J — 1), the elemental mass and stiffness matrices are 



Me = 





/ 


156 


22h 


54 


-13h \ 


h 




22h 


Ah^ 


13h 




420 




54 


13h 


156 


-22h 




I 


-13h 


-3/i2 


~22h 


4/i2 J 



2P 
1^ 





6 


3h 


-6 




3/1 


2h^ 


-3h 




-6 




6 


V 


3/1 




-3/i 



3/i 

-3h 



and the global system H4.1|l is obtained as an assembly of the previous elemental 
matrices. 

For the resolution of 1)4. l|l at each time step, we use the following lemma with 
A={M + At^pS) ,X^At^ and / = F" 

Lemma 4.1. 0, [HI Let A be a symmetric positive definite 2 J x 2 J real matrix, 
f G M^'^ and u' be the solution of Au' = f . Then, for all A > 0, the system 



with Kh 



p2,7~2 



Au + \d^^{u) 3 /, 
X [~g, g] X R admits an unique solution u given by 

{Au)^ =/,;, ^r^-l,...,2J-2,2J. 



We consider a steel pipe of length L = 1.501 m, with an external diameter equal 

EI 

to 1 cm and a thickness equal to 0.5 mm. Thus — — = 282.84 m'^.s~^ where 

pS 



E — 2x 10^^ Pa is the Young's modulus, p = 8 x lO'^ kg/m^ is the material density, 
S is the cross-section and / the cross-sectional moment of inertia of the pipe (see 
also [7] for a more detailed description of the mechanical setting). The vibration 
of the support is given by (f){t) = 0.2sin(10t) for alH > 0, 17 = 0.1 and the initial 
data are uq = 0, vq = —2h i.e u(-, 0) = ut{-, 0) = (at < = the beam is at rest). 

In the next figures we show the approximate motion of the impacting end of the 
beam. The results given at figures El and ^ have been obtained with /3 = 1/2, 
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J = 19 and At = 5 X lO'^s, At = lO^^s and At = 5 x lO^^s (let us recall that we 
have unconditional stability for this value of 



fem-Newmark scheme 




0.0 0.2 0.4 0.6 0.8 1.0 i.2 1.4 i.6 1.8 2.0 



Figure 2. At = 5x IQ-^s and /3 = i 




We can observe that the trajectories are almost the same at the beginning of 
the time interval (up to the end of the first "contact period" i.e < t < 0.2) and 
remain quite similar afterwards even if the details of the impact phenomenon are 
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Figure 4. At = IQ-^s and /3 = i 



different. This is not surprising since vibrations with unilateral constraints always 
lead to sensitivity to initial data. 

The motion of the impacting end of the beam has also been computed by using 
the normal compliance approximation of Signorini's conditions. In this case we 
have to define the penalty parameter e. Although the corresponding stiffness 1/e 
has a physical meaning, the range of values usually chosen is quite large: 1/e — 
lO^^A^.m"^ in 9 , 1/e = 5.5 x lO^iV.rn"^ in ^3] for instance. In the following 
results we consider 1/e = lO^N.m^^ and we apply once again a Newmark's scheme 
with J = 19 and Ai = 5 X 10~S, At = 10'^ s and At = 5 x IQ-'^s (see figures 
|5||7I). We should notice that we have to solve now a partial differential equation, 
thus we choose /3 = 1/4 for which the unconditional stability of Newmark's scheme 
holds. 



Nevertheless we observe a kind of numerical instability (figure EJ: spurious high 
frequencies appear during " contact periods" and this phenomenon can be controlled 
only for very small time steps. Moreover, the non-penetration condition is violated 
by the trajectories computed with the normal compliance approximation while it 
remains satisfied by the approximate motions uf^ ^. More precisely we have 

max \uavv(L, t) — n\ = 2.25 x 10^**m for At = 5 x 10^'^s 

0<t<2' ' 

which is a rather coarse approximation of the unilateral constraint. 

Finally let us point out that, for J = 19 and (3 = 1/4, the stability condition 
(prT|l leads to 



At < 2. -^ = 3.3469 x lO^S 
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Normal Compliance - Newmark scheme 




0.000 0.194 0.387 0.581 0.775 Om9 1. 162 1.356 1.550 



Figure 5. At = 5 x 10~S and /3 = i 



Normal Compliance - Newmark scheme 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Figure 6. At = lO'^s and (3 = I 

with k{h) given by 1)5. 2|l (see the Appendix). Thus, with /? = 1/4 we can solve 
the discretized problems iPhi3^)o<n<N-i with J — 19 and At = 5 x 10~^s (fig- 
ure ISJ: the computation remains less expensive than with the normal compliance 
approximation (see the next table for a comparison of the CPU times). 



CPU 
Time 



Fig. 2 Fig. 3 Fig. 4 Fig. 5 Fig. 6 Fig. 7 Fig. 8 



18.24 173.71 870.45 241.17 1187.9 2426.4 859.38 
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Normal Compliance - Nevnnark schei 




0.0 0.2 U.4 0.6 O.S 1.0 1.2 1.4 1.6 1.8 2.0 



Figure 7. At = 5 x IQ-'^s and /3 = i 



fern- trapezoidal , 




Figure 8. At = IQ-^s and /3 = i 



Remark 4.2. The computations has been performed on MAC G4 (1.1 Ghz) with 
Scilab, the scientific computing software developed by ENPC-INRIA. Other results 
are available at the following web address: ,http:/ /www.univ-reunion.fr~ydumont/beam-vibrations.htm( 



5. Appendix: estimate of K{h). 

We consider the finite element space discretization described at section 01 i.e 
Vh = span{iy9i, . . . ,f2j} where {ipi)i<i<2j is the Hermite piecewise cubics basis. 
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Thus, for all Uh G Vh we have 

.7 J 
Uh = ^Uh (x^) ip2i-i + Uft (a;,) ^Pn- 

i=l i=l 

In order to simplify the notations, we let 

Ui = Uh {xi) u'i = u'f^ {xi) for alH = 0, . . . , J. 

We may observe that, since Uh & Vh C V , we have uq = u'q = 0. 
Let us compute first {uh,Uh)- We have 



13 . , 9 , 



i=i 



Then we use the following relations 



^ X Wj-i'^o-i - u'jUj) Aa; = - ^ {Axf u'jUj, 



and we get 



K, ^Z') = ^ E (y ("j- + + A^ K-1 - "j)) A^ + 7^ E + "i-i) A^ 

i=i i=i 

1 / A / ^ n2 ^ 1752 ^ \ 294 2 A 

Then, we observe that 

1 2 294 

— (Ax«:;-9^.,;) A:r + ^z.,Aa: 

= I T^Axuj - ^uj\ Ax + -^UjAx + fJi^'^L {u'jf Ax 
420 V 19 2 J 3024 420 x 19^ ^ 
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and thus 

,/ 



140 

1T52 ^ — > / 2 2 \ A 1T52 .9 X — ^ / \2 A 

+7560 g + + rmo^""-'-' + 180 g - 

2 



1 /18^ , 19 51 2^ 

— Axu'j uj Ax H uiAx. 

420 V19 2 y 3024 



We deduce the following inequality 
(5.1) 



(Ax)^ 5^ ((^)' + K-i)') Aa; + ^ K - Ax 



> 



420 X 192 , 
Now we compute a {uh, Uh)- 



J J 



a {uh,Uh) = 4 ( 12 E ("j + "j-i) Aa; - 24 Y^MjUj_iAx 
(Ax) \ 



+ ^4 I 12Ax^ (uj-iu'j — Uju'j_i) Ax + 4 (Aa;)^ ^ Wj-u^.^Aa; 
I 4 (Ax)' ^ ((«;.)' + (4_0') Ax - 12 (Ax)' uju'j I , 



(Ax) 
which gives 



J 



v2 



(Ax')EK -«;-i) Ax|. 



Hence 



.7 



a{uh,uh) < I 24^(uj-Uj_i)'Ax + 6(Ax2)^(^+^_i)'Ax 



[Axf 

(Ax^)EK-^-i)'Ax 



0=1 j=i 
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Using H5.1|l . we deduce 

24fc2 420 X 192 

a{uh,Uh) < 4 — iuh,uh), VuheVh 

(Ax) -J' 

and thus 

, , 24 X 420 X P 2^x5x192x21 P 

(5.2) K ih) < — 



37 {Axf 37 (Ax)^' 

Remark 5.1. This is certainly not an optimal upper bound for k (h). 
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